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Abstract. We calculate decadal aerosol direct and indi- 
rect (warm cloud) radiative forcings from US anthropogenic 
sources over the 1950-2050 period. Past and future aerosol 
distributions are constructed using GEOS-Chem and histor- 
ical emission inventories and future projections from the 
IPCC A1B scenario. Aerosol simulations are evaluated 
with observed spatial distributions and 1980-2010 trends of 
aerosol concentrations and wet deposition in the contiguous 
US. Direct and indirect radiative forcing is calculated us- 
ing the GISS general circulation model and monthly mean 
aerosol distributions from GEOS-Chem. The radiative forc- 
ing from US anthropogenic aerosols is strongly localized 
over the eastern US. We find that its magnitude peaked in 
1970-1990, with values over the eastern US (east of 100° W) 
of — 2.0 Wm -2 for direct forcing including contributions 
from sulfate (— 2.0 Wm -2 ), nitrate (— 0.2 Wm -2 ), organic 
carbon (—0.2 W m -2 ), and black carbon (+0.4 W m -2 ). The 
uncertainties in radiative forcing due to aerosol radiative 
properties are estimated to be about 50 %. The aerosol in- 
direct effect is estimated to be of comparable magnitude to 
the direct forcing. We find that the magnitude of the forc- 


ing declined sharply from 1990 to 2010 (by 0.8 Wm -2 di- 
rect and l.OWm -2 indirect), mainly reflecting decreases in 
SO 2 emissions, and project that it will continue declining 
post-2010 but at a much slower rate since US SO 2 emis- 
sions have already declined by almost 60 % from their peak. 
This suggests that much of the warming effect of reducing 
US anthropogenic aerosol sources has already been realized. 
The small positive radiative forcing from US BC emissions 
(+0.3 W m -2 over the eastern US in 2010; 5 % of the global 
forcing from anthropogenic BC emissions worldwide) sug- 
gests that a US emission control strategy focused on BC 
would have only limited climate benefit. 


1 Introduction 

Growth in population and energy demand over the past 

100 yr has greatly increased the anthropogenic source of at- 

mospheric aerosols in the United States. This has caused 

public health, visibility, and deposition concerns (US En- 

vironmental Protection Agency (US EPA), 2009, 2010). 


Published by Copernicus Publications on behalf of the European Geosciences Union. 


3334 


E. M. Leibensperger et al.: Climate response 


Increasingly strict regulations on aerosol sources have been 
enacted by US air quality agencies over the past decades. 
However, aerosols also exert a negative radiative forc- 
ing on climate and reductions in their abundance aggra- 
vate greenhouse-driven climate change (Raes and Seinfeld, 

2009) . We quantify in this paper the changes in radiative 
forcing arising from historical and projected trends of US 
anthropogenic aerosol sources for the 1950-2050 period. In 
a companion paper (Leibensperger et al., 2012), we use a 
general circulation model (GCM) to analyze the resulting cli- 
mate response. 

Anthropogenic aerosols mainly consist of sulfate, nitrate, 
ammonium, black carbon (BC), and organic carbon (OC). 
Sulfate and nitrate aerosols are formed by oxidation of SO 2 
and nitrogen oxides (NO x = NO + NO 2 ). Coal combustion is 
the dominant source of SO 2 in the US. SO 2 emissions grew 
until 1980 and then decreased by 56% between 1980 and 
2008 (US EPA, 2010). Anthropogenic NO x is emitted by 
fuel combustion in general and US emissions decreased by 
36% between 1990 and 2008 (US EPA, 2010). Ammonia 
originates mainly from agriculture (Bouwman et al., 1997) 
and has not been subjected to regulation. BC is emitted from 
small, low-temperature combustion sources such as residen- 
tial, transport, and small industrial fuel burning. Organic 
aerosol is traditionally partitioned by models into a primary 
component (POA) from combustion and a secondary compo- 
nent (SOA) of dominant biogenic origin (Kanakidou et al., 
2005). According to Bond et al. (2007), US anthropogenic 
emissions of BC and POA decreased from 1925 to 1970 due 
to a decline in residential coal use, but increased from 1970 
to 1990 due to increases in broader fuel use. 

Aerosols directly affect climate by scattering and absorb- 
ing solar radiation. The Intergovernmental Panel on Cli- 
mate Change (IPCC) estimates the aerosol direct effect to 
presently exert a global mean negative radiative forcing 
of —0.5 zb 0.4 Wm -2 , partly offsetting the positive radia- 
tive forcing of +2.6 zb 0.3 W m -2 from the long-lived green- 
house gases (Forster et al., 2007). Aerosols also indirectly 
affect climate by modifying cloud properties. Absorbing 
aerosols affect cloud cover by increasing solar heating in 
the atmosphere (“semi-direct effect”; Koch and Del Genio, 

2010) . Aerosols can also act as cloud condensation nuclei 
(CCN) and ice nuclei. An increase in CCN reduces cloud 
droplet size, which brightens the clouds (“cloud albedo ef- 
fect”; Twomey, 1974) and enhances cloud cover by reducing 
precipitation efficiency (“cloud lifetime effect”; Albrecht, 
1989). The IPCC best estimate for the global indirect radia- 
tive forcing from the cloud albedo effect is —0.7 W m -2 , with 
a range of uncertainty from —0.3 to — 1.8 Wm -2 (Forster 
et al., 2007). The cloud lifetime effect is even more un- 
certain, but could be comparable in magnitude to the cloud 
albedo effect with additional effects on the hydrological cy- 
cle (Lohmann and Feichter, 2005; Denman et al., 2007). Ad- 
ditional aerosol indirect radiative effects involving cloud ab- 


sorption, height, or glaciation could also be significant (Den- 
man et al., 2007). 

The need to integrate air quality and climate change mit- 
igation objectives in environmental policymaking is increas- 
ingly recognized (National Research Council, 2005; Raes 
and Seinfeld, 2009; Penner et al., 2010). This is particularly 
the case for aerosols since air quality improvements poten- 
tially come at the cost of warming. Previous studies have in- 
vestigated the climate effects of global aerosol sources (Shin- 
dell et al., 2008; Chen et al., 2010b; Koch et al., 2011), but 
understanding the effects of national sources is most useful 
for policymakers. The US is an interesting testbed to analyze 
the climate implications of environmental regulations since 
the historical period from 1950 to present has witnessed a 
reversal of aerosol trends, increasing until 1980 and then de- 
creasing, with regulations in place to enforce continued de- 
crease in the future. Recent GCM simulations by Mickley 
et al. (2012) suggest that completely removing US anthro- 
pogenic aerosol sources would increase temperatures in the 
eastern US by 0.4-0. 6 °C on an annual mean basis and as 
much as by 1-2 °C during summer heat waves. 

Unlike other aerosol components, BC exerts a positive di- 
rect radiative forcing on climate. It has been argued that de- 
creasing BC emissions (and hence aerosol absorption) could 
provide a “win-win” strategy for air quality and climate 
change mitigation (Jacobson, 2002; Bond, 2007; Grieshop 
et al., 2009). However, such a strategy is complicated by the 
fact that BC sources cannot be controlled in isolation from 
other aerosol components, in particular POA (Unger et al., 
2010). There is additional uncertainty regarding the sign of 
the semi-direct and indirect effects associated with BC (Koch 
and Del Genio, 2010; Koch et al., 2011; Chen et al., 2010a). 

We use here a global chemical transport model (GEOS- 
Chem CTM) to reconstruct historical aerosol trends from 
1950 to present and project future trends to 2050, with fo- 
cus on the US aerosol loadings and evaluation with observed 
trends. We then use these monthly mean aerosol distributions 
as input to the Goddard Institute for Space Studies (GISS) 
GCM 3 to construct a 1950-2050 timeline of aerosol direct 
and indirect radiative forcing with decadal resolution, resolv- 
ing the contributions from the different aerosol components. 
The aerosol distributions presented here are used in a com- 
panion paper (Leibensperger et al., 2012) to determine the 
sensitivity of US climate to anthropogenic aerosol sources in 
the past and in the future. 

2 Methods 

2.1 Global aerosol simulation 

We use GEOS-Chem CTM simulations of coupled tro- 
pospheric ozone-NO x -VOC-aerosol chemistry (version 
8.01.01; http://geos-chem.org/) to describe the global evo- 
lution of aerosol concentrations resulting from changes in 
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anthropogenic emissions from 1950 to 2050. The simula- 
tions are conducted for a series of 2 -yr decadal time slices 
from 1950 to 2050. The first year is used for initialization 
and the second year for analysis. All simulations use the 
same 2000-2001 meteorological data from the NASA 
Goddard Earth Observing System (GEOS-4). The GEOS-4 
data are available with 1 x 1.25° horizontal resolution, 
55 levels in the vertical, and a temporal resolution of 6 h 
(3 h for surface variables). The data are regridded here to 
2 x 2.5° horizontal resolution for input to GEOS-Chem. 
Using the same meteorological year for the 1950-2050 
simulations isolates the effects of emission changes. Climate 
change over that period might affect monthly mean aerosol 
concentrations by ± 0 . 1-1 pgm -3 , with even the sign of 
the effect uncertain (Jacob and Winner, 2009). Aerosol 
abundances are more sensitive to the large emission changes 
simulated here. 

GEOS-Chem simulates the aerosol mass concentrations of 
sulfate-nitrate-ammonium (SNA), POA, SOA, BC, and size- 
resolved sea salt and soil dust (Park et al., 2006; Liao et al., 
2007). Sea salt and soil dust are viewed as natural and are not 
discussed here further. Gas-phase and aerosol chemistry are 
coupled by in-cloud SO 2 oxidation, gas-aerosol thermody- 
namic partitioning of SNA and SOA, aerosol effects on pho- 
tolysis rates (Martin et al., 2003), and heterogeneous chem- 
istry (Jacob, 2000; Evans and Jacob, 2005). Water-soluble 
gases and aerosols are scavenged in convective updrafts as 
well as by rainout and washout from convective anvils and 
large-scale precipitation (Liu et al., 2001). The model allows 
for species to return to the atmosphere if falling precipitation 
evaporates. Dry deposition is modeled using a resistance in 
series scheme (Wesely, 1989) as implemented by Wang et 
al. (1998). 

Sulfate is formed by gas-phase reaction of SO 2 with OH 
and by aqueous-phase reaction of SO 2 with H 2 O 2 and ozone 
in clouds and sea-salt aerosols (Alexander et al., 2005). The 
global mean tropospheric lifetime of sulfate in the model 
(computed as the ratio of global burden to deposition) is 
4.0 days, comparable to other sulfate models (Schulz et al., 
2006). Formation of SNA aerosol is computed locally with 
the MARS-A aerosol thermodynamic equilibrium model 
(Binkowski and Roselle, 2003). Following Park et al. (2003), 
we assume that 20 % of BC and 50 % of POA emitted from 
anthropogenic sources is hydrophilic and thus available for 
in-cloud scavenging. The remaining portions of BC and POA 
are emitted as hydrophobic and become hydrophilic in the 
atmosphere with an e-folding time of 1.2 days. Formation 
of SOA in GEOS-Chem involves gas-aerosol partitioning of 
semi-volatile oxidation products from isoprene and terpenes 
(Heald et al., 2006; Liao et al., 2007). Global production of 
SOA is 32 Tg C a -1 , similar in magnitude to previous GEOS- 
Chem studies (Henze et al., 2008; Pye et al., 2010). 

Figure 1 shows the 1950-2050 trends of global and US 
aerosol sources. Emissions of SO 2 and NO x are from 
EDGAR Hyde 1.3 (van Aardenne et al., 2001) for 1950- 



1950 2000 2050 1950 2000 2050 

Year Year 


Fig. 1. Global (solid) and US (dashed) trends in emissions of SO 2 , 
NO x , BC, and POA for 1950-2050. US emissions are multiplied by 
10 to fit on scale. SO 2 and NO x emissions are from EDGAR (van 
Aardenne et al., 2001; Olivier and Berdowski, 2001). BC and POA 
emissions are from Bond et al. (2007). All emissions are extended 
past the year 2000 following the IPCC A1B scenario. 


1990 and EDGAR 3.2 FT (Olivier and Berdowski, 2001) for 
the year 2000. Emissions of BC and POA are from Bond 
et al. (2007). 2010-2050 decadal emissions are calculated 
by applying growth factors to year 2000 emissions. Similar 
to Fiore et al. (2002) and Wu et al. (2008), growth factors 
for different categories of anthropogenic emissions (biofuel 
and fossil fuel) and countries are derived from the Integrated 
Model to Assess the Greenhouse Effect (IMAGE; Streets et 
al., 2004) following the IPCC A1B scenario (Nakicenovic 
and Swart, 2000). The decline of US anthropogenic SO 2 
and BC in the A1B scenario is comparable to those in the 
more recent RCP6.0 scenario from the IPCC AR5 (Moss et 
al., 2010). All RCP scenarios decrease global SO 2 , NO x , 
BC, and POA and US NO x emissions more rapidly than the 
A IB scenario. We do not consider trends in global emis- 
sions from aircraft (0.07 TgS a -1 ; Chin et al., 2000) and 
ships (4.2 Tg S a -1 ; Corbett et al., 1999). 

Natural sources of SO 2 , NO x , and VOCs are held constant 
for the 1950-2050 period. Volcanic emissions of SO 2 (non- 
eruptive and continuously eruptive only) are from Andres 
and Kasgnoc (1998) and amount to 5.5TgSa -1 . Oceanic 
dimethylsulfide (DMS) in GEOS-Chem has a global source 
of 21 TgSa _1 and produces 19TgSa -1 of SO 2 following 
oxidation by OH and NO 3 (Park et al., 2004). Lightning 
NO x emissions are calculated from a parameterization based 
on cloud top height (Price and Rind, 1 992) and then scaled 
locally to match satellite observations from the Optical Tran- 
sient Detector/Lightning Imaging Sensor (Sauvage et al., 
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2007; Murray et al., 2012). The global emission of light- 
ning NO x is 5.6TgNa _1 . Soil NO x emissions follow the 
algorithm of Yienger and Levy (1995) as implemented by 
Wang et al. (1998); this amounts globally to 6.6TgNa _1 
including 0.7TgNa _1 from fertilizer application (also held 
constant over 1950-2050). Biogenic emissions of isoprene 
and monoterpenes are calculated using the Model of Emis- 
sions of Gases and Aerosols from Nature (MEGAN; Guen- 
ther et al., 2006), which yields globally 380TgCa _1 iso- 
prene and 98TgCa _1 monoterpenes. We use the climato- 
logical biomass burning inventory from Duncan et al. (2003), 
which yields l^TgSa -1 of S0 2 , 6.5TgNa _1 of NO x , 
2.9 TgC a -1 of BC, and 23 Tg C a -1 ofPOA. 

Ammonia emissions are also held constant over the 1950- 
2050 period for lack of better information. We use the 
global inventory of Bouwman et al. (1997) overwritten in 
Asia by Streets et al. (2003). Seasonal variations of am- 
monia emissions are described by Park et al. (2004) and 
are constrained in the US to match wet deposition data 
(Gilliland et al., 2003). Global emissions total 60TgNa _1 
including 40TgNa -1 from anthropogenic sources (mainly 
agriculture), 14TgNa -1 from natural biogenic activity, and 
5.9TgNa -1 from biomass burning. US emissions total 

2.2 TgN a -1 from anthropogenic sources, 0.6 TgN a -1 from 
natural biogenic activity, and 0.04 TgN a -1 from biomass 
burning. 

2.2 Aerosol direct radiative forcing 

We use the GISS GCM 3 (Rind et al., 2007) as modified by 
Chen et al. (2010b) to calculate all-sky aerosol direct and 
indirect radiative forcing. The GCM uses monthly mean tro- 
pospheric aerosol distributions from GEOS-Chem including 
sulfate, nitrate, BC, and OC (OC = POA + SOA). OC is mul- 
tiplied by 1 .4 to convert from the simulated mass of organic 
carbon to total organic matter (Malm et al., 1994). Aerosol 
water content is calculated using hygroscopic growth fac- 
tors from the Global Aerosol Data Set (GADS; Koepke et 
al., 1997). These aerosol components are incorporated in a 
climate equilibrium simulation with resolution of 4° latitude 
x 5° longitude, 23 vertical levels extending from the surface 
to 0.002 hPa, and fixed present-day sea surface temperatures 
(SST) and sea ice (Rayner et al., 2003). All-sky aerosol di- 
rect radiative forcing is determined within the GCM through 
parallel radiative calculations including and excluding an- 
thropogenic aerosols. “Anthropogenic” includes contribu- 
tions from fuel use and industry, but not open biomass burn- 
ing. Sea salt and soil dust concentrations are from Hansen 
et al. (2002) and do not influence the direct radiative forcing 
calculated here. 

Radiative forcing calculations are conducted for both ex- 
ternal and internal aerosol mixtures. Externally mixed 
aerosol components are assumed to have a standard gamma 
size distribution with an area-weighted effective variance of 
0.2 and an effective dry radius of 0.3 pm for sulfate and ni- 


trate, 0.5 pm for OC, and 0. 1 pm for BC (Chung and Seinfeld, 
2002; Liao et al., 2004). Refractive indices are from Toon 
et al. (1976) for sulfate and nitrate, and from d’ Almeida et 
al. (1991) for BC and OC. As in Chung and Seinfeld (2002) 
and Liao et al. (2004), we determine the aerosol extinction 
efficiency, single scattering albedo, and asymmetry parame- 
ter for each aerosol component from a lookup table generated 
by offline Mie calculations as a function of aerosol size and 
refractive index. Radiative forcing is calculated by perturb- 
ing each aerosol component individually. The total anthro- 
pogenic aerosol radiative forcing is calculated by perturbing 
all aerosol components. 

Internally mixed aerosols are composed of sulfate, nitrate, 
OC, and BC. As in Chung and Seinfeld (2002), we assume 
that the water uptake of internally mixed particles is equal 
to the sum of water uptake for externally mixed particles. 
The internally mixed particles are assumed to have a stan- 
dard gamma size distribution with surface area- weighted dry 
radius of 0.3 pm and area-weighted variance of 0.2 (Chung 
and Seinfeld, 2002; Liao et al., 2004; Chen et al., 2007). 
Optical properties of the internal mixture are calculated us- 
ing the volume-weighted mean of the refractive indices of 
the individual components. The radiative forcing of an indi- 
vidual aerosol component is calculated by removing it while 
holding the remaining components fixed. 

2.3 Aerosol indirect radiative forcing 

Chen et al. (2010b) previously applied the GISS GCM 3 to 
calculate the aerosol indirect effects initiated by aerosol par- 
ticles acting as CCN and thus altering the number concentra- 
tion N c of cloud droplets. We follow their approach here us- 
ing the gridded monthly aerosol concentration fields of SNA 
and carbonaceous aerosols from GEOS-Chem, together with 
their sea salt aerosol concentrations. N c is related to the con- 
centration of water soluble aerosol ions, mi, by a standard 
power law dependence (Boucher and Lohmann, 1995; Pen- 
ner et al., 2006): 

log A c = A + B logmi (1) 

We calculate N c from the archived GEOS-Chem aerosol dis- 
tributions. Gridded 3-D monthly mean fields of A and B 
were obtained by Chen et al. (2010b) from detailed simula- 
tions of sulfate and sea salt aerosol microphysics and acti- 
vation within the GCM (Adams and Seinfeld, 2002; Nenes 
and Seinfeld, 2003; Fountoukis and Nenes, 2005; Pierce 
and Adams, 2006). While the microphysical simulations fo- 
cused on sulfate and sea salt aerosols, the relationship be- 
tween the coefficients (A and B) and N c is general and ad- 
ditionally applicable to nitrate and hydrophilic carbonaceous 
aerosols. Chen et al. (2010b) applied these coefficients to 
global aerosol fields and we use the same A and B fields 
here for each simulation. Small changes to A and B may 
be expected as aerosol composition evolves, but we do not 
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address this issue here. As in Chen et al. (2010b), N c is re- 
stricted to be 20 cm -3 or greater. This restriction reduces 
radiative forcing estimates in pristine environments (Chen et 
al., 2010b) but has little effect in the US, our area of interest. 

Following Chen et al. (2010b), we consider the cloud 
albedo and lifetime effects applied to liquid stratiform clouds 
only. The cloud albedo effect arises from the enhancement 
of cloud optical depth from the presence of smaller cloud 
droplets. The cloud optical depth scales as the inverse of the 
area-weighted mean effective radius r e of the cloud droplet 
size distribution (Del Genio et al., 1996). r e is related to N c 
by 


where L is the liquid water content of the cloud (cm 3 wa- 
ter per cm 3 air), and k is a constant (0.67 over land, 0.80 
over ocean, Martin et al., 1994) that relates the volume mean 
droplet radius and r e . We restrict r e to be greater than 2 pm 
and less than 20 pm. 

The cloud lifetime effect involves the effect of N c on 
the rate of autoconversion of cloud droplets to precipitation. 
Chen et al. (2010b) introduced a dependence of the autocon- 
version rate on N c based on the work of Khairoutdinov and 
Kogan (2000), which fitted results from large eddy simula- 
tions of drizzling stratocumulus: 

|! = -1350 rq^N-'™ (3) 

where q\ is the cloudwater mass content (kg water per kg of 
air) and y is a tuning parameter added by Hoose et al. (2008) 
and Chen et al. (2010b) to the original Khairoutdinov and 
Kogan (2000) equation (y = 1) in order to retain GCM cli- 
mate equilibrium. We find climate equilibrium is retained 
for y = 12. This value is consistent with the factor of 2 to 10 
underestimate of autoconversion rates found to occur when 
using gridbox-scale values of N c (Morales andNenes, 2010). 

The radiative forcing from the cloud albedo effect is cal- 
culated for each decade between 1950 and 2050 using par- 
allel radiative calculations, similar to the calculation of the 
aerosol direct radiative forcing. The parallel radiative cal- 
culations are conducted using cloud optical properties de- 
termined from N c distributions with and without US an- 
thropogenic aerosols. In this method, the cloud albedo ef- 
fect is calculated without the complications of climate feed- 
backs and can be treated as a traditional radiative forcing. 
The cloud lifetime effect cannot be calculated in this man- 
ner due to coupling with the hydrological cycle. Instead, 
multi-year climate equilibrium simulations are required to 
account for the effects of feedback processes on cloud distri- 
butions (Denman et al., 2007; Forster et al., 2007). We thus 
calculate the total aerosol indirect radiative forcing (cloud 
albedo and lifetime effects) for each decade between 1950 


and 2050 by conducting two 60-yr climate equilibrium sim- 
ulations, one with and one without US anthropogenic aerosol 
sources. Both simulations use fixed greenhouse gases, SST, 
and sea ice for year 2000. The first 10 yr are used for initial- 
ization and the following 50 yr (representing climate equi- 
librium conditions) are used for analysis. Comparison to the 
radiative forcing of the cloud albedo effect, which is calcu- 
lated independently as described above, allows separate as- 
sessment of the first and second aerosol indirect effects. In 
order to test the significance of our results, we perform a Stu- 
dent’s t-test accounting for autocorrelation as described by 
Zwiers and von Storch (1995). 

Our use of monthly mean aerosol distributions archived 
from GEOS-Chem does not allow for feedbacks between 
aerosols and cloud properties (Chapman et al., 2009). These 
feedbacks are likely very small relative to the changes of 
aerosol sources studied here. Using monthly mean aerosol 
distributions does not significantly bias our estimates of the 
direct radiative effect (Koch et al., 1999), but likely influ- 
ences our calculation of the aerosol indirect effects due to 
the nonlinear relationship between aerosol amount and cloud 
droplet number (Jones et al., 2001). 


3 Evaluation of 1980-2010 US aerosol trends 

We evaluate our simulation of aerosol sources, concentra- 
tions, and related long-term trends over the US with 1980- 
2009 wet deposition flux data available from the National 
Atmospheric Deposition Program/National Trend Network 
(NADP/NTN; http://nadp.sws.uiuc.edu/NTN), and 1990- 
2009 surface air concentrations available from the US Inter- 
agency Monitoring of Protected Visual Environments (IM- 
PROVE; http://vista.cira.colostate.edu/IMPROVE) and the 
Clean Air Status and Trends Network (CASTNET; http: 
//] ava.epa.gov/castnet). The wet deposition flux data provide 
constraints on aerosol sources considering that most of SO 2 , 
NO x , and ammonia emitted in the US is deposited within the 
country by wet and dry processes (Chin and Jacob, 1996; Li 
et al., 2004; Zhang et al., 2012). 

Figure 2 compares the simulated and observed annual wet 
deposition fluxes and surface concentrations of sulfate, ni- 
trate, and ammonium across the US for the years 1980, 1990, 
2000, and 2010. Observations are three-year average values 
centered on the decadal year except 2010, which is a two- 
year average for 2008-2009. All observation sites meeting 
the NADP/NTN data completion criterion for at least two of 
the three years are included. Observations of surface air con- 
centrations are not continuous and some records have sub- 
stantial gaps. To address this, we first aggregate the obser- 
vations as seasonal means (DJF, MAM, JJA, SON) for each 
year, requiring for each season a minimum of 10 data points 
for IMPROVE (3 -day averages) and 5 data points for CAST- 
NET (weekly averages). The seasonal means for the three 
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Fig. 2. Sulfate wet deposition fluxes (kg ha -1 a -1 , left) and sul- 
fate concentrations in surface air (pgm -3 , right). Observations 
(circles) are 3-year averages for 1979-1981 (wet deposition only), 
1989-1991, 1999-2001, and 2008-2009. GEO S-Chem model val- 
ues (background contours) are from decadal time slice simulations 
with 1980, 1990, 2000, and 2010 emissions. Reduced major axis 
linear regressions between model and observations are shown inset 
for each map along with the coefficient of determination (r 2 ) and 
normalized mean bias (NMB). 


years are then averaged, with the additional requirement that 
data be available for at least two of the years. 

We assemble comparison statistics between the model (M) 
and observations ( O ) using the coefficient of determination 
(r 2 ), the corresponding reduced-major-axis linear regression, 
and the normalized mean bias (NMB) calculated for all N 
sites containing valid data: 


NMB = 


E,=l (Ml - Oj) 
ZL o, 


( 4 ) 


We additionally calculate the mean normalized bias MNB = 
E,=iKV - Oi)/Oil The MNB is generally within 5 % 
of the NMB and smaller than the NMB. 

Figure 2 shows that the model has excellent agreement 
with observed sulfate wet deposition across the US, captur- 
ing the spatial variability (r 2 = 0.71-0.78) and magnitude. 
The NMB is less than 4 % for all available decades except 
1980 (—15%). Simulated surface concentrations of sul- 


fate have similar success, capturing spatial variability (r 2 = 
0.79-0.95) with low bias (10 % or less). 

Figure 3 shows the simulated and observed trends in sul- 
fate wet deposition and air concentrations over the 1980- 
2010 period for the eastern US (east of 100° W). Values 
are averaged over all sites ( n , shown inset) with observa- 
tions available for more than 25 yr (deposition) and more 
than 17 yr (concentrations). Observed wet deposition of sul- 
fate decreased by 58% between 1980 and 2010, consistent 
with the 56% reduction in SO 2 emissions reported by US 
EPA (2010). Sulfate aerosol concentrations show a paral- 
lel decreasing trend for 1990-2010. Modeled deposition and 
surface concentrations both decrease by 40% over 1980- 
2010, less than observed because of an apparent underesti- 
mate of 1980 emissions. The simulated decreases of sulfate 
deposition and surface concentrations are in better agreement 
with observations after 1990. 

Figure 4 compares model and observations for ammonium 
wet deposition fluxes and aerosol concentrations. The model 
successfully captures the distribution (r 2 = 0.63-0.75) and 
magnitude (NMB < 7 %) of ammonium deposition between 
1980 and 2010. The largest deposition is in the agricultural 
Midwest where emissions are highest. In that region the 
model is lower than observed. Aerosol ammonium concen- 
trations peak in the industrial Midwest, reflecting the combi- 
nation of high ammonia emissions and high concentrations 
of sulfate and nitrate that promote the fractionation of am- 
monia into the aerosol as ammonium. The model captures 
this spatial distribution (r 2 = 0.72-0.84), but is biased high, 
particularly in 2010 (+42 %). 

Figure 3 shows observed and simulated 1980-2010 trends 
in ammonium wet deposition and surface concentrations in 
the eastern US. The wet deposition data show no significant 
trend, consistent with the model (which assumes constant 
ammonia emissions). However, the observed surface con- 
centrations show a large decreasing trend that is only weakly 
reproduced in the model. This decrease can be explained by 
decreases in sulfate and nitrate under conditions when the 
aerosol is neutralized. The inability of the model to capture 
the decrease in ammonium concentrations reflects a problem 
with the nitrate trend, as discussed below. 

Figure 5 shows the wet deposition and surface air concen- 
trations of nitrate. The model captures the spatial distribution 
of nitrate deposition (r 2 = 0.71-0.81), but underpredicts to- 
tal nitrate wet deposition by 20-33 % in 1980-2000. The low 
bias vanishes in 2010, following a decade of large reductions 
of NO x from the energy generation sector (US EPA, 2010). 
Aerosol nitrate in the model is exclusively ammonium ni- 
trate, and its formation is contingent on ammonia availability 
(since formation of ammonium sulfate is favored thermody- 
namically over ammonium nitrate). As such, its simulation 
tends to compound errors in sulfate and ammonium (Park 
et al., 2006). Nevertheless, our simulation shows relatively 
little bias and is improved compared to previous versions 
of GEOS-Chem (Park et al., 2004, 2006). The observed 
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Fig. 3. 1980-2010 trends in wet deposition fluxes (solid) and surface aerosol concentrations (dashed) in the eastern US (east of 100° W). 
Values are annual means. Model trends (blue) are compared to observations (black). The number of sites averaged in each calculated trend 
is given inset. Model values are sampled at the site locations and then averaged. 
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Fig. 4. Same as Fig. 2 but for ammonium. 


Fig. 5. Same as Fig. 2 but for nitrate. 


maximum over the Midwest is well captured by the model 
where it reflects a balance between ammonia, NO x , and SO 2 
source influences. 

Observations of nitrate wet deposition show a 33 % de- 
crease from 1980 to 2009, consistent with the reported 36 % 


reduction of NO x emissions (US EPA, 2010). Most of this 
reduction was realized between 2000 and 2009 through the 
NO x State Implementation Plan (SIP) Call. However, the 
model is largely flat over 1980-2010. The EDGAR emis- 
sions used for 1980-2000 are apparently too low. The 
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IMAGE A1B projection used for 2000-2010 does not ac- 
count for the NO x SIP Call, and this fortuitously produces 
a good match to the 2008-2009 observations of nitrate wet 
deposition (Fig. 5). 

Observed aerosol nitrate concentrations were flat from 
1990 to 2000 but then decreased by 23 % between 2000 and 
2009. In contrast, model concentrations increased weakly 
during this period because of the decrease in sulfate, result- 
ing in more ammonium being available for ammonium ni- 
trate formation (West et al., 1999). The observed 2000-2009 
decrease in nitrate, rapid decrease in sulfate, and near con- 
stant supply of ammonia seems at odds with current under- 
standing of SNA thermodynamics. We did not explore this 
issue further, but it warrants future study. 

Figure 6 compares simulated and observed surface con- 
centrations of BC and OC for 1990-2010. Biogenic VOCs 
are a major source of OC in the Southeast. Open fires are an 
important component of observed variability for BC and OC 
in the West and in the Southeast (Park et al., 2007). These 
fires are not well reproduced by the model on an event ba- 
sis, explaining in part the relatively low correlations between 
model and observed concentrations. The model is also biased 
low by 15-40% on average. Previous GEOS-Chem studies 
did not show as severe a bias because the BC and OC sources 
were fitted to reproduce the observations (Park et al., 2003) 
and included significantly higher biomass burning emissions 
(Park et al., 2003, 2006). Previous models using the Bond et 
al. (2004) inventory (similar to Bond et al. (2007) used here) 
also underestimate BC in North America (Koch et al., 2007, 
2009). 

Figure 3 shows that observed BC and OC concentrations 
over the US decreased by 50% and 34% respectively be- 
tween 1990 and 2009. The model trends are much weaker, 
with BC and OC concentrations decreasing by 27% and 
16% respectively, as driven by the trends in the Bond et 
al. (2007) emission inventory and the IMAGE A IB projec- 
tion. The mismatch between simulated and observed trends 
reduces the model low bias in 2010 to 20 % for BC and 15 % 
for OC. 

Overall, the model shows mixed success in the simulation 
of 1990-2010 aerosol trends. We have an excellent simula- 
tion for sulfate, which accounts for 65 % of anthropogenic 
aerosol mass and 78% of associated aerosol optical depth 
over the eastern US (Table 1). This gives us some confidence 
in the model representation of overall trends in aerosol ra- 
diative forcing, and we will quantify below the related error 
using the trend information in Fig. 3. In view of our un- 
derestimate of BC for the 1990-2000 period, we will focus 
our discussion of BC radiative forcing on present-day (2010) 
when the simulation shows no significant bias. 


Black Carbon Organic Carbon 



0.00 0.25 0.50 0.75 1.00 0.00 0.75 1.50 2.25 3.00gg rrr 3 

Fig. 6. Black carbon (left) and organic carbon (right) concentration 
in surface air (pgm -3 ). Observations (circles) are 3-yr averages 
for 1989-1991, 1999-2001, and 2008-2009. GEOS-Chem model 
values (background contours) are from the decadal time slice simu- 
lations with 1990, 2000, and 2010 emissions. Reduced major axis 
linear regressions between model and observations are shown inset 
for each map along with the coefficient of determination (r 2 ) and 
normalized mean bias (NMB). 


4 Aerosol direct radiative forcing from US 
anthropogenic sources 

Figure 7 (top) shows the calculated annual mean aerosol 
direct radiative forcing for the peak in US anthropogenic 
aerosols in year 1980. The values in Fig. 7 assume an internal 
aerosol mixture, the forcing increases by 10% for an exter- 
nal mixture due to the lesser effect of BC absorption (Chylek 
et al., 1995; Jacobson, 2000). The forcing is strongly lo- 
calized over the eastern US, where it reaches — 4.9Wm -2 . 
This more than offsets for that region the 1980 radiative 
forcing from the long-lived greenhouse gases (+1.8 W m -2 ). 
On a global scale the aerosol direct radiative forcing from 
US anthropogenic aerosol sources in 1980 is — O.lOWm -2 , 
25 % of our computed global aerosol direct radiative forcing 
of —0.41 Wm -2 from anthropogenic sources worldwide in 
1980. By 2010 we find that the radiative forcing from anthro- 
pogenic US aerosol sources has decreased to —0.03 W m -2 
globally, amounting to just 6% of the total from world- 
wide anthropogenic sources (— 0.54 Wm -2 ), reflecting the 
rapid decline of emissions in the US and growth in Asia 
(Fig. 1). The decline of the US contribution to the global 
mean aerosol radiative forcing results largely from decreas- 
ing sulfate aerosol (+0.07 W m -2 ) with a smaller opposing 
role from decreasing black carbon (—0.01 W m -2 ). 
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Table 1. Aerosol direct radiative forcing from global and US anthropogenic sources (2010) a . 



Emissions (Tg a 1 

Burden (Gg) c 

AOD x 1000 d 

Direct Radiative Forcing 
External Mixture (W m -2 ) 

Direct Radiative Forcing 
Internal Mixture (W m -2 ) 

Global 

| Global 

Eastern US 

Global 

Eastern US 

Global 

Eastern US 

Global 

Eastern US 

Sulfate 

80.3 (8%) 

390 (6.5 %) 

8.9 (65%) 

34.2 (6.6 %) 

83.2 (67%) 

-0.65 (6.9 %) 

-1.85 (66%) 

-0.62 (7.0%) 

-1.82 (64%) 

Nitrate 

35.6(15%) 

46.1 (2.0%) 

1.1 (58%) 

6.1 (2.7%) 

18.1 (65%) 

-0.11 (2.7%) 

-0.35 (63 %) 

-0.10(2.9%) 

-0.33 (63 %) 

Black Carbon 

4.3 (7.1%) 

69.1(6.3%) 

1.5 (66%) 

0.9 (6.6%) 

1.9 (64%) 

+0.19(5.5%) 

+0.32 (65 %) 

+0.30 (5.0%) 

+0.49 (67 %) 

Primary Organic Aerosol 

8.1 (5.5%) 

110(6.1%) 

2.3 (69 %) 

1.4 (6.7%) 

3.2 (69%) 

-0.14(6.1%) 

-0.25 (67 %) 

-0.11 (5.5%) 

-0.20 (68 %) 

Total 

- 


- 

42.6(6.1%) 

106(66%) | 

-0.67(6.5%) 

-2.13(65%) 

-0.54(6.1%) 

-1.88 (64%) 


a Model values for anthropogenic perturbations (excluding biomass burning) averaged globally and over the eastern US (east of 100° W). The contribution from US anthropogenic 
sources is shown in parentheses. 

b Sulfate emissions as precursor SO 2 (Tg S a -1 ); nitrate emissions as precursor NO* (Tg N a -1 ); black carbon (Tg C a -1 ); primary organic aerosol (Tg C a -1 ). 
c Sulfate burden as Gg S; nitrate burden as Gg N; black carbon burden as Gg C; organic carbon burden as Gg C. 
d Aerosol optical depth (AOD) computed for an external mixture. 


Table 1 presents the 2010 anthropogenic perturbations to 
aerosol abundances and optical properties averaged over the 
globe and over the eastern US (east of 100° W). We see 
that the US contributes only 5 % of the global mean anthro- 
pogenic radiative forcing from absorbing aerosols (global 
mean: +0.30 Wm -2 ) while providing 6.3% of the forcing 
from scattering aerosols (global mean: — 0.83 Wm -2 ). We 
show in Leibensperger et al. (2012) that even though the 
aerosol radiative forcing from US anthropogenic sources is 
very small on a global scale, the large regional forcing elicits 
a strong regional climate response over the eastern US. 

Figure 8 (top) shows the 1950-2050 evolution of the 
annual mean aerosol direct radiative forcing from US an- 
thropogenic sources averaged over the eastern US (east of 
100° W). It peaks in 1970-1990 at — 2.0 Wm -2 (internal 
mixture). When externally mixed, the radiative forcing is 
stronger, — 2.2 Wm -2 . The external mixture has 5% lower 
column integrated aerosol optical depth, 2% higher single 
scattering albedo, and 1 % higher asymmetry parameter (for- 
ward scattering is increased) over the eastern US. In the ex- 
ternal mixture the radiative forcing of BC is 40 % lower and 
that of sulfate and nitrate is 5-7 % higher. OC shares com- 
mon sources with BC, which makes its contribution to radia- 
tive forcing more sensitive to mixing state; it has 30 % higher 
radiative forcing in an external mixture than an internal mix- 
ture. 

The anthropogenic radiative forcing trends in Fig. 8 mirror 
the model US emission trends in Fig. 1 . We previously evalu- 
ated the model trends for 1980-2010 using observed records 
of wet deposition and aerosol concentrations (Fig. 3). For 
sulfate, which is the dominant anthropogenic radiative forc- 
ing agent, the model closely reproduces the observed trends 
for 1990-2010, but is 15% too low in 1980, suggesting a 
corresponding error in the aerosol radiative forcing estimate 
for 1980. For BC, the model is 40% too low in 1990 and 
2000 but the bias decreases to 20% in 2010, reflecting a 
steep decline in the observed concentrations. Radiative forc- 
ing from OC is relatively small due to its lower extinction ef- 
ficiency and modest anthropogenic component. The anthro- 
pogenic component of OC simulated here consists mainly of 
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Fig. 7. Annual mean aerosol direct (top), first indirect (middle), 
and total indirect (bottom) radiative forcing from US anthropogenic 
sources for year 1980. The aerosol direct effect is calculated assum- 
ing an internal mixture. Assuming an external mixture would in- 
crease the radiative forcing over the US by 10 %. The direct and first 
indirect effects are calculated with parallel radiative calculations us- 
ing perturbed aerosol abundances (direct) and cloud optical prop- 
erties (first indirect). The total aerosol indirect effect is calculated 
by difference between 50-year climate equilibrium simulations with 
and without US anthropogenic aerosols. White areas show no sig- 
nificant change at the 95th percentile of confidence. The global 
mean values are —0.07 W m -2 for the direct effect, —0.05 W m -2 
for the first indirect effect, and — 0.08 W m -2 for the total indirect 
effect. 

POA with limited changes in SOA. Radiative forcing from 
nitrate is weaker than its annual mean concentration would 
imply because it is most abundant in winter when insolation 
is weak. 
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Aerosol Direct Forcing from US Anthropogenic Aerosols 



Year 

Aerosol Indirect Radiative Forcing from 
US Anthropogenic Aerosols 


0.5 



- 2.5 - 


- 3.0 

1 950 1 960 1 970 1 980 1 990 2000 2010 2020 2030 2040 2050 
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Fig. 8. 1950-2050 evolution of aerosol direct (top) and total (first + 
second) indirect radiative forcing (bottom) from US anthropogenic 
sources over the eastern US (east of 100° W). Circles represent the 
net direct radiative forcing and the bars give the contributions from 
the different components. The direct radiative forcing calculations 
are for an internal aerosol mixture. Aerosol indirect radiative forc- 
ing is calculated as the difference between two 50-yr equilibrium 
climate simulations with vs. without US anthropogenic aerosols. 
The error bars represent the 95 % confidence interval of the forcing 
based on model interannual variability. 


We can use the data in Fig. 3 to estimate the model error 
in computing the change in anthropogenic aerosol direct ra- 
diative forcing from 1980 to 2010. Our best model estimate 
for this change (adding the contribution from the different 
aerosol components, black line in Fig. 8) is +0.78 Wm -2 . 
Applying scaling factors for errors in the individual compo- 
nents based on Fig. 3, and assuming that the model bias for 
BC in 1990 applies also to 1980, we obtain a change in an- 
thropogenic radiative forcing constrained by observations of 
+0.87 W m -2 , about 10 % larger. In addition to model biases 
in aerosol loading, uncertainties surround the aerosol radia- 
tive properties required to connect aerosol concentrations to 
radiative forcing. Schulz et al. (2006) compared the radiative 


forcing estimates of nine models using identical emission 
datasets. The relative error in radiative forcing attributable 
to aerosol optical properties varied slightly between aerosol 
components but was about 50 %. We adopt this value as the 
uncertainty of our radiative forcing calculations. 

Future projection of US emissions from the IMAGE model 
applied to the IPCC A1B scenario (Fig. 1) indicates a contin- 
ued decrease in aerosol direct radiative forcing (Fig. 8), but 
with little change past 2020. Aerosol radiative forcing will 
have decreased by almost a factor of 4 from its peak strength 
in 1970-1990. 

Two policy-relevant implications can be drawn from the 
trends in US aerosol direct radiative forcing in Fig. 8. First, 
future decrease in aerosol radiative forcing due to reduc- 
tions in anthropogenic emissions will likely be smaller than 
the decrease that has already been realized over 1980-2010. 
Thus most of the climate response from controlling US an- 
thropogenic aerosol sources has already been realized, and 
almost all will have been realized by 2020. Second, the 
present-day (2010) radiative forcing from US anthropogenic 
BC is small (and even less if external forcing is assumed), so 
that emission controls targeting BC provide only limited cli- 
mate benefit. These emission controls would benefit public 
health but are not an effective “win-win” scenario for both 
US public health and climate. The climate benefit of reduc- 
ing US BC emissions is further complicated when consider- 
ing the effects of co-emitted aerosol species and cloud inter- 
actions (Bauer et al., 2010; Chen et al., 2010a; Koch and Del 
Genio, 2010; Unger et al., 2010). We elaborate further on 
the implications of US anthropogenic BC in Leibensperger 
et al. (2012). 

5 Aerosol indirect radiative forcing from US 
anthropogenic sources 

Estimating the aerosol indirect radiative forcing from cloud 
albedo and lifetime effects is far more uncertain than the di- 
rect radiative forcing. Uncertainties relate to the conversion 
from aerosol mass concentration to cloud droplet number 
concentration (N c ) and from there to cloud optical proper- 
ties, including the effective cloud droplet radius (r e ) and the 
spatial distribution of liquid water path (vertical column of 
liquid water amount, gm -2 ). Particularly uncertain is the 
subgrid variability of these parameters. Table 2 compares 
the global mean cloud properties computed from our model 
to Chen et al. (2010b) (same GCM and cloud droplet param- 
eterization, but different aerosol concentrations) and to other 
values in the literature. Our values are consistent with the 
published ranges. Chen et al. (2010b) report higher global 
mean N c but similar cloud forcing, which may reflect differ- 
ences in the geographical and vertical distributions of N c . 

Figure 9 shows the change in liquid stratiform cloud prop- 
erties due to US anthropogenic aerosols in 1980 including 
column N c , effective cloud droplet radius, and cloud optical 
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Table 2. Global mean cloud properties* 1 . 



This 

Work 

Chen et 
al. (2010b) 

Other 

Studies b 

Total Column N c (10 10 m 2 ) 

3.8 

6.2 

2.1-7. 6 C 

N c at 850 hPa (cm“ 3 ) 

83.0 

122.2 

75. 0-135. 0 d 

r e at cloud top (pm) e 

12.9 

12.2 

6.8— 13.3 f 

LWP (g m -2 ) 

111.3 

111.5 

41.5-1 10. 0 g 

TOA SW CF (W m -2 ) 

-55.2 

-53.8 

-46.0— 61. 0 h 


a Cloud droplet number concentration (N c ), effective area- weighted mean cloud 
droplet radius (r e ), liquid water path (LWP), top-of-atmosphere shortwave cloud 
forcing (TOA SW CF). 

b Ranges from the literature including models and observation analyses. 
c Menon et al. (2002), Gettelman et al. (2008), Hoose et al. (2008), and Salzmann et 
al. (2010). 

d Penner et al. (2006). 

e Modeled values of r e are multiplied by 2+ 3 to approximate the cloud top conditions 
of satellite retrievals (Meskhidze et al., 2007). 

f Menon et al. (2002), Kristjansson et al. (2005), Penner et al. (2006), Storelvmo et 
al. (2006), Gettelman et al. (2008), and Barahona et al. (2011). 
g Kristjansson et al. (2005), Penner et al. (2006), Storelvmo et al. (2006), Gettelman 
et al. (2008), Hoose et al. (2008), and Salzmann et al. (2010). 
b Menon et al. (2002), Kristjansson et al. (2005), Penner et al. (2006), Gettelman et 
al. (2008), and Salzmann et al. (2010). 


depth. US anthropogenic aerosols increase column inte- 
grated N c by up to 8 x 10 10 m -2 (40%) which reduces the 
effective cloud droplet radius of liquid stratiform clouds by 
up to 1.2 pm (15%), and increases the warm cloud optical 
depth by more than 3 (17%). The resulting radiative forc- 
ing from the cloud albedo effect amounts to 67% of the 
direct forcing (— 1.3 Wm -2 averaged over the US east of 
100° W) on an annual basis (Fig. 7 middle) and reaches up 
to — 3.0 Wm -2 in the Southeast. Similar to the aerosol di- 
rect effect, radiative forcing from the cloud albedo effect is 
largely confined to the US and the downwind North Atlantic. 

Figure 7 (bottom) shows the total aerosol indirect radiative 
forcing (cloud albedo and lifetime) from US anthropogenic 
aerosols in 1980. The indirect forcing is calculated as the 
mean difference in net top of the atmosphere radiation be- 
tween 50-yr simulations with and without US anthropogenic 
aerosol sources, as described in Sect. 2.3. It is noisier than 
the direct radiative forcing, largely reflecting the cloud life- 
time effect, but is similarly concentrated over the eastern 
US and the North Atlantic downwind. The annual mean 
total indirect forcing in the eastern US (east of 100° W) is 
—2.2 W m -2 for 1980, 10 % larger than the direct effect and 
70 % larger than the cloud albedo effect alone. The total indi- 
rect effect of US anthropogenic aerosols is also slightly larger 
globally (— 0.08 Wm -2 total indirect vs. — 0.07 Wm -2 di- 
rect). 

Figure 8 (bottom) shows our computed 1950-2050 trend 
in aerosol indirect radiative forcing over the eastern US (east 
of 100° W) from US anthropogenic sources. The error bars 
indicate the 95 % confidence interval of this difference based 
on model interannual variability. We estimate this uncer- 
tainty to be zb0.5Wm -2 . The trend in indirect forcing is 
consistent with the evolution of SO 2 emissions (Fig. 1) and 
aerosol direct forcing (Fig. 8 top). Indirect forcing peaks in 
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Fig. 9. Effect of US anthropogenic aerosols on column integrated 
cloud droplet number concentration, liquid stratiform cloud droplet 
effective radius at cloud top (r e ), and liquid stratiform cloud optical 
depth. Values are annual model means for 1980. 

1970-1990 at a value of — 2.0 Wm -2 , the same magnitude 
as the aerosol direct effect. SO 2 regulations cause a sharp 
decrease in forcing with a change of +1.0 Wm -2 between 
1990 and 2010, and little change afterward. 

6 Conclusions 

Efforts to improve air quality through aerosol source reduc- 
tions could have significant regional climate implications 
due to the strong and localized radiative forcing exerted by 
aerosols. The US is of particular interest for investigating this 
effect as aerosol concentrations peaked in the 1970-1990 pe- 
riod, have decreased rapidly since then, and are projected to 
continue decreasing in the future. We used here the GEOS- 
Chem CTM applied to historical emission inventories and 
future projections (IPCC A IB scenario) to simulate global 
trends in aerosol concentrations over the 1950-2050 period, 
and we applied those in the GISS GCM to calculate the di- 
rect and indirect radiative forcings from US anthropogenic 
aerosol sources. A companion paper by Leibensperger et 
al. (2012) examines the resulting climate response over the 
US and globally. 

A prominent feature of the aerosol radiative forcing his- 
tory from US sources is the rapid 1980-2010 decline driven 
by air quality regulations. We tested our reconstruction of 
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this history with observed spatial distributions and long-term 
trends in wet deposition (1980-2010) and aerosol concentra- 
tions (1990-2010). We reproduce successfully the observed 
45 % decline of sulfate wet deposition and concentrations be- 
tween 1990 and 2010, but underestimate the earlier 1980 to 
1990 decline. Ammonium wet deposition fluxes show no 
long-term trend in the observations or in the model (where 
ammonia emissions are assumed constant). The combination 
of decreasing SO 2 emissions and constant ammonia emis- 
sions leads to a weak increasing trend of nitrate in the model, 
but observations show a decrease. Observations of BC and 
OC show stronger 1980-2010 declines than simulated by the 
model, implying a faster decline of emissions than is repre- 
sented in current inventories. 

We calculated a global direct radiative forcing from US an- 
thropogenic aerosol sources of —0.10 Wm -2 in 1980 when 
the sources were at their peak. This forcing is strongly lo- 
calized over the eastern US and downwind North Atlantic. It 
averages — 2.0 Wm -2 over the eastern US (east of 100° W), 
including — 2.0 Wm -2 from sulfate, +0.4 Wm -2 from BC, 
—0.2 W m -2 from nitrate, and —0.2 W m -2 from OC. These 
values assume an internal aerosol mixture; an external mix- 
ture would increase the overall negative radiative forcing by 
10% due to decreased absorption by BC. Based on prior 
modeling studies, our radiative forcing estimates have an un- 
certainty of about 50 % due to the uncertainty of aerosol ra- 
diative properties. 

The direct radiative forcing from US anthropogenic 
aerosols declined strongly between 1970-1990 and 2010, 
largely reflecting the decline in SO 2 emissions. The global 
forcing declined to -0.03 W m -2 by 2010, amounting to 
only 6 % of the forcing from global anthropogenic aerosol 
sources (as compared to 25 % in 1980). The forcing over the 
eastern US declined by 0.78 Wm -2 in the model between 
1990 and 2010. Correction to this model estimate based on 
the observed trends of wet deposition fluxes and concentra- 
tions would imply a decline of 0.87 W m -2 . 

We calculated the aerosol indirect forcing including the 
first and second indirect effects (cloud albedo and cloud life- 
time) applied to warm stratiform clouds. US anthropogenic 
aerosols potentially alter properties of convective and ice 
clouds but the scientific understanding of these effects is very 
low (Denman et al., 2007). We found that the indirect radia- 
tive forcing from US anthropogenic aerosol sources is mainly 
localized over the eastern US and North Atlantic, similarly to 
the direct forcing, and shows similar temporal trends over the 
1950-2050 period. While uncertainties are large, 1980 val- 
ues over the eastern US average — 2.2 Wm -2 . The indirect 
forcing is similar in magnitude to the direct radiative forcing. 
About 60 % of the indirect radiative forcing comes from the 
cloud albedo effect. 

We project that by 2050 the direct and indirect radiative 
forcings over the eastern US from US anthropogenic sources 
will have decreased to — 0.6 Wm -2 and — 0.9 Wm -2 , re- 
spectively, 71 % and 55 % weaker than the 1980 peak, and 


46 % and 25 % weaker than present-day (2010). Nitrate is ex- 
pected to become a major component of the aerosol radiative 
forcing in 2050 as NO x emissions decrease more slowly than 
SO 2 emissions. Essentially all of the 1980-2050 decrease in 
radiative forcing occurs over the 1990-2020 period. US an- 
thropogenic emissions are sufficiently weak by now that little 
additional reduction in forcing is expected in the future even 
with the projected continued decrease in emissions. 

The positive radiative forcing from BC anthropogenic 
emissions in the US is presently too small to provide sig- 
nificant climate leverage from future emission controls tar- 
geted at BC. Our simulations significantly underestimate US 
BC concentrations in earlier periods but we have confidence 
in our estimate of the present-day radiative forcing from US 
anthropogenic BC because of the ability of the model to re- 
produce observed BC concentrations in the US in 2010 (20 % 
low bias). Although our radiative forcing estimate does not 
include the semi-direct effects associated with cloud evap- 
oration (Koch and Del Genio, 2010), transient-climate sim- 
ulations presented in our companion paper (Leibensperger 
et al., 2012) do include these effects and confirm the cli- 
mate insensitivity to present-day anthropogenic sources of 
BC in the US. This insensitivity merely reflects the relatively 
small magnitude of BC emissions from the US and does not 
challenge the argument that decreasing global BC emissions 
would have significant benefit for both air quality and climate 
(Jacobson, 2002; Bond, 2007; Grieshop et al., 2009; Shindell 
et al., 2012). In particular, an important regional climate ben- 
efit could be achieved in Asia, where BC sources are much 
larger than in the US (Lu et al., 201 1). 
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